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Abstract. We analyze the structure of the X-ray emis- 
sion of a sample of 22 Abell clusters of galaxies with a cD 
in their centre, observed with the ROSAT PSPC. Utilizing 
the multi-scale power of the Wavelet Transform we detect 
significant (~ 50 h" 1 kpc) offsets between the large scale 
centroid and the peak of X-ray emission. Despite the un- 
certainties on the satellite pointing, the X-ray to optical 
correlation indicates a likely association between the X- 
ray peak and the dominant galaxy. We develop a model 
in which the offset is produced by small amplitude oscil- 
lations of the cD galaxy around the bottom of the cluster 
potential well and successfully compare it to the observed 
distribution of offsets. Within this scenario the offsets are 
not due to dynamic instabilities and the number of struc- 
tured clusters is greatly reduced. 
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1. Introduction 

It is well known that the presence of a dominant galaxy 
(cD) in the core of a galaxy cluster is a sign of dynamical 
evolution, being highly correlated with a smooth spatial 
distribution and a low spiral fraction of member galaxies 
(Sarazin 1988). 

In a large fraction of clusters with cDs, the X-ray peak 
is offset with respect to the centroid of large scale emis- 
sion (Mohr et al. 1993, 1995). Whether these offsets are 
caused by substructure, due to the infall of smaller groups 
of galaxies, or not, is still an open issue. Indeed, the frac- 
tion of clusters with substructure is strongly linked to the 
density parameter fig (Richstone et al. 1992; Nakamura 
et al. 1995) and any spurious substructure would confuse 
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the potential link between observed cluster morphologies 
and that expected theoretically. 

In a sample of 22 X-ray images of Abell clusters with 
cD galaxies in their centre (see section || for the descrip- 
tion of the sample), we find that the peak of the X-ray 
emission is always located about the position of the dom- 
inant galaxy and is offset from the large scale centroid. In 
this work we propose an alternative to infall or Intra Clus- 
ter Medium (hereafter ICM) instability models to explain 
the presence of these offsets since: a) these phenomena 
are unexpected in clusters which have the appearance of 
being relaxed systems (cooling flows and dominant galax- 
ies); b) infall would disrupt the cooling flow and produce 
substructure with higher frequency in the outskirts of clus- 
ters rather than in their centre, and c) the asymmetries 
observed in the ICM would have to be independent of the 
cD position. 

In our model a strong and compact X-ray source is 
associated with the dominant galaxy which in turn moves 
about the bottom of the cluster potential well with har- 
monic oscillations, producing the observed offset. Whether 
this compact source is a cooling flow or an active nucleus 
is not essential and will be discussed in the last section. 
Such a model provides a natural explanation of the facts 
summarized above without invoking any dynamical insta- 
bility. As a consequence, the fraction of structured clus- 
ters is reduced from more than 70% (Davis 1994) down to 
about 30%. 

The role of cD galaxies and their motion relative to 
the bottom of the cluster potential well is still a matter of 
debate. Optical observations, which should be best suited 
for such an analysis, suffer from low number statistics in 
the cluster and the results are somewhat contradictory. 
Studies on individual clusters (Hill et al. 1992; Sharpies et 
al. 1988; Oegerle & Hill 1992) find significant peculiar ra- 
dial motions of the cD with respect to the mean recession 
velocity, and in sample of 25 cD clusters Zabludoff et al. 
(1990) find that "a substantial fraction of the cD galax- 
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Table 1. Properties of the cluster sample and of the respective PSPC observations. 
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800518 
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800249 
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94.0 
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9794 
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8881 


369.0 






OUUOZD 








2657 


0.0414 


800320 


18904 


12487 


44.0 


3562 


0.0490 


800237 


20199 


13503 


0.0 


3921 


0.0936 


800378 
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0.0292 


800354 


3353 (3) 


4845 




4059 


0.0488 


800175 
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' ' Approximate number of photons (background subtracted) used in the large scale centroid determination 
' 2 ' Cooling flow mass rate in solar masses per year (from Fabian 1994; White et al. 1997; Allen & Fabian 1998). A long dash 
indicates clusters for which no data have been found. 

' 3 ' This exposure is lower than the 5 ks limit set for the sample, but a sufficiently accurate measurement of the centroid has 
been possible. 



ies have velocities significantly different from the mean of 
their parent clusters" , a result confirmed by Malumuth et 
al. (1992). However it is not completely clear whether the 
presence of substructure in the velocity field could alter 
this result (Oegerle & Hill 1994; Bird 1994). 

Numerical simulations of cluster evolution (Malumuth 
& Richstone 1984) predict that cD galaxies are created 
and spend their existence near the centre of clusters and 
move slowly around them. Moreover, the actual motion of 
the dominant galaxy and the related displacement from 
the bottom of the potential well are also an important 
clue to understand the formation scenario of dominant 
galaxies (see Oegerle & Hill 1994 and references therein). 

By assuming that cD galaxies move about the centre 
of mass of their parent clusters we can derive an equation 
of motion and estimate a normalized oscillation amplitude 
which is independent from the parameters of the individ- 
ual clusters. We can then compare the expected distri- 
bution of the normalized amplitudes to the measurements 
from the set of 22 X-ray images of moderate redshift Abell 
clusters. 

This paper is organized as follows: section b| describes 
the data and the reduction procedures, section yand || de- 
scribe respectively the small and large scale analysis per- 
formed on the images, while section ^| briefly summarize 
the measured offset properties. The possible identification 



of the cD galaxy with the emission peak is discussed in sec- 
tion [| section describes the physical oscillatory model 
and the results arc summarized and discussed in section pi 



2. The data: selection criteria and preliminary re- 
duction 

To detect a small but significant offset between the peak 
and the large scale emission of clusters we need a precise 
determination of both these positions. 

The uncertainty on the determination of the position 
of a source characterized by a profile with a width a (irre- 
spective of whether it is intrinsic or due to the instrument 
resolution) and a number N of counts is given by: 



A(a:) ~ -j= 
y/N 



(1) 



which follows from the standard error propagation rules. 
The exact equality holds for a Gaussian spatial distribu- 
tion of photon counts. 

Hence, for a point-like source, the best instrument to 
measure the position maximizes the value of the ratio 
yfe/cTfjy where e is the detector efficiency and its 
intrinsic resolution. However, when the angular dimension 
of the source is intrinsic - as is the case for the ICM X-ray 
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emission - the width of the photon distribution is fixed and 
only a higher value of e can increase the position accuracy. 

Given the above discussion, the cluster observations 
have been extracted from the ROSAT/PSPC archive 
maintained by the ROSAT group at the Max Plank In- 
stitut fur Extraterrestrische Physik in Garching (D). We 
have selected a set of observations pointed on Abell clus- 
ters with cD galaxy in their centre and a known value of 
the radial velocity dispersion. Other constraints were an 
exposure time greater than 5 ks and an angular dimension 
an order of magnitude greater than the resolution in the 
centre of the PSPC field (~ 40" FWHM). 

As of March 1995, we found 22 clusters with 
ROSAT/PSPC pointed observations that met the above 
requirements, listed in table [j] along with their main prop- 
erties. All these observations have JV^5x 10 3 counts in 
the cluster large scale emission. With a typical width of 
3r c ~ 7.5' for the clusters as a whole, equation [l] gives an 
accuracy of 6" in the determination of the centroid, where 
r c is the core radius of the modified King function that 
better fits the cluster profile. We have used 3r c as cluster 
width to properly take into account the wider wings of the 
King with respect to the Gaussian distribution. 

The central compact source - which in our model is 
responsible for the emission peak - has a mean number 
of ^ 1000 photons and a width comparable to the PSPC 
central resolution of ~ 25" (la). Equation [j] gives an ac- 
curacy A(x) £ 1". 

The precision on the determination of the peak posi- 
tion could be increased using a detector with higher reso- 
lution as the ROSAT HRI. However, the efficiency of this 
instrument is lower than that of the PSPC and hence the 
accuracy of centroid measurement would be decreased. 
Note that to avoid any influence of the boresight errors 
(see section ^3j) it is crucial that both the centroid and the 
peak positions are measured on the same image. 

The image reduction has been performed with the 
complete treatment included in the ESAS software (Snow- 
den 1995) which takes into account all the properties of the 
PSPC instrument and of the particular observation (solar 
contamination and short and long term enhancements of 
the cosmic background). This is because instrumental ef- 
fects such as vignetting, detector efficiency, particle back- 
ground etc. conspire to build a large scale structure in the 
images. This structure, if uncorrected, has its centre coin- 
cident with the image centre and could efficiently mimic 
an offset with the peak of the true emission. 

We have tested the complete reduction process search- 
ing for large scale residuals in PSPC point-like source im- 
ages. We found no evidence for such an effect in the anal- 
ysis of a set of 20 images. 

Due to different distance and intrinsic cluster dimen- 
sions, we have implemented specific ESAS routines in or- 
der to obtain images with dynamic pixel size. The fi- 
nal images (512 x 512 pixels each) have a pixel size (1 



xVd.o.f. = 1.16 : 




Centroid shift (pixel) 

Fig. 1. Influence of the compact source on the centroid deter- 
mination. Histogram of the distance between input and output 
centroids are plotted for the simulations with (shaded region) 
and without (solid line) the offset compact source. The distri- 
butions are consistent with being two realizations of the same 
parent. See the text for the statistical test results. 

pixel = 4 -j- 10") chosen to maintain a fixed (small) ratio 
between the whole image width and the cluster angular 
dimensions. 

To fully explore any possible bias due to the reduction 
process we could simulate a set of cluster images with 
instrumental effects included and than compare the clus- 
ter parameters obtained after the reduction with those 
input in the simulations. However, this procedure is ex- 
tremely time consuming and any unknown instrumental 
effect would remain undetected. 

3. Small scale analysis 

Once a cleaned count-rate image has been obtained, a 
wavelet transform algorithm is applied to search for small 
scale structures embedded in the cluster emission. We use 
the algorithm expressly designed to detect and charac- 
terize small structures overlaid on strongly varying back- 
grounds presented and discussed in Lazzati et al. (1998). 

The wavelet transform algorithm spans scales from half 
the PSF on-axis up to about twice the PSF at the edges 
of the image, with a spacing of a factor of two between ad- 
jacent scales. This procedure guesses the position of the 
emission peak and selects a set of "point sources" to be 
masked in the subsequent large scale analysis (see sec- 
tion |). 

The emission peak position is then refined fitting a 
bidimensional Gaussian source model to the image in the 
wavelet space. Since the covariance matrix highly under- 
estimates uncertainties in correlated data sets, the uncer- 
tainties on the peak position are estimated in x 2 space. A 
X 2 shift of 4.61 is used in order to obtain 90% confidence 
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Fig. 2. Influence of the compact source on the centroid de- 
termination. Measured x-axis errors on the centroid position 
versus the photon counts in the offset peak are plotted in the 
upper panel. The vast majority of measurement have a disper- 
sion lower than one pixel, but a limited number of centroid 
measurement are shifted towards the compact source (positive 
values in both panels) . No clear dependence of this shift on the 
power of the compact source is seen. The lower panel show the 
same position errors in the measured centroid versus the offset 
of the input compact source. In this case a trend with the com- 
pact source centroid is visible. The trend disappears for offsets 
larger than 20 pixels, almost the size of the core radius of the 
large scale emission. Note that in both cases the errors in the 
centroid measurements tend to mask the effect of displacement 
between the centroid and the emission peak. 



intervals for two independent degrees of freedom (position 
in both directions). Again, this procedure is not strictly 
correct in the presence of correlations between data points. 
Rosati (1995) performed several tests to see if any sys- 
tematic error could arise, confirming the reliability of this 
procedure. 

Note, however, that the uncertainties on the peak po- 
sition are in most cases negligible with respect to those 
in the large scale centroid, which contribute the largest 
uncertainties to the offset measurement. 



4. Large scale analysis 

Great care must be used in the determination of the cen- 
troid position since the offset we are looking for is consid- 
erably smaller than the whole cluster emission. To obtain 
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KS = 22% \ 
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Fig. 3. Histogram of the measured errors in recovering the 
large scale centroid normalized to the uncertainties derived 
from the covariance matrix. The distribution appears consis- 
tent with a Gaussian distribution with a normal standard de- 
viation (shaded area). 

the best position for the centroid of the large scale emis- 
sion, we fit the cluster emission with an elliptical (3 model: 



Ix (r, i?) oc 



1 



-3/3+| 



(2) 



where r c (i?) is the azimuth dependent core-radius. The 
central part, affected by the presence of the compact sour- 
ce, was excluded from the fit as well as all the point-like 
sources detected in projection over the cluster emission. A 
bidimensional Gaussian smoothing of the image was nec- 
essary to obtain a more accurate estimate of the statistical 
distribution of points. 

To test the reliability of this procedure we have sim- 
ulated a set of 1000 clusters with properties that span 
the range of our sample which are summarized in ta- 
ble | The core radius of the clusters has been fixed to 
20 pixelsEl while all the remaining shape parameters have 
been randomly picked up each time from the distributions 
described in table ||. To fully reproduce the experimental 
images, an offset Gaussian source [FWHM = 12 pix- 
els, slightly bigger than a point-like source) with a lower 
number of counts has been added in the cluster centre 
and then the centroid has been measured with the pro- 
cedure described above. Analogous simulations have been 
performed without the compact offset Gaussian source for 
comparison. 

Figure [j], || and || show the result of these simulations. 
To see whether the presence of the compact offset source 
could produce a systematic error in the centroid deter- 
mination we have compared the position errors measured 



1 This is not a restriction of simulations since the image size 
has been chosen to preserve a fixed dimension of the cluster in 
pixels. 
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Table 2. Properties of the clusters simulated to test the cen- 
troid determination procedure. All distributions are flat except 
ellipticity and /3 which have Gaussian probability functions. 



LS (1) 



(2) 



f3 SS (3) A(» (4) A(y) 



(4) 



average 15000 0.8 0.76 5000 40 40 
width 10000 0.07 0.25 3000 20 20 
^' Number of counts in the large scale cluster emission. 

Ellipticity defined as the ratio between the minor and 
major axis. 

Number of photon counts in the offset compact source. 
^ Offset of the compact source in the x and y directions in 
pixels with respect to the large scale centroid. 



from simulations with and without the compact source. 
Figure [l] shows the histogram of the distance between in- 
put and output centroid positions in both cases. The sim- 
ilarity between the two distributions is confirmed by the 
Kolmogorov-Smirnov test, which gives a ~ 21% proba- 
bility of them both deriving from the same parent dis- 
tribution. Moreover, by fitting the difference of the two 
distributions with a zero constant we have obtained a re- 
duced x 2 value of 1.16. The absolute values of the distance 
between the true and the measured centroid are around a 
quarter of a pixel, which, given the resolution of our im- 
ages, translates in an uncertainty in the range 1" — 2.5", 
a value suitable for our requirements (see section ^|). 

When the offset Gaussian source is present in simula- 
tions, the distribution of the shift between input and out- 
put centroids has a tail that is considerably higher than 
what expected in a Gaussian distribution (the positive 
values in figure f|. This effect increases with the shift be- 
tween the centroid and the peak position (lower panel of 
figure ||) up to a sharp cutoff when the shift becomes larger 
then the core radius. We can confidently ignore this effect 
since it tends to reduce the measured offset and not to 
create it spuriously. 

Finally, we have tested the reliability of the uncertain- 
ties extracted from the covariance matrix. If these are cor- 
rect, the errors measured in the simulations divided by 
those estimated in the fitting procedure should be dis- 
tributed as a Gaussian with normal standard deviation. 
As shown in figure || this is what actually happens. The 
observed distribution is consistent with a Gaussian parent 
at the 22% level with a Kolmogorov-Smirnov test. A fit 
with a Gaussian gives a reduced \ 2 of 1.17. 



5. Offsets 

Applying the procedure described and tested above to the 
X-ray images we are able to measure the offsets of the 
cluster sample along with their uncertainties. The values 
of centroid position errors estimated in the fitting proce- 
dure were always 5" but, to be conservative, we adopt 
this value for the whole set of clusters. 
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Fig. 4. Positions of the large scale centroids with respect to 
the emission peak. The peak is located in the origin of coordi- 
nates and shifts in right ascension and declination are given in 
arcseconds. The dashed circle encloses the 5" confidence region 
(see text) for the centroid position while dotted circle marks 
the 10" region. 



The results are shown in figure [| where the emission 
peak lie at the center of coordinates and the symbols mark 
the position of the centroid. The inner circle radius is 5" 
(1 a), while the outer is 10". Only 4 measurements over 22 
are in agreement with the absence of any offset at the 2 a 
level. Figure || shows the offsets distribution in physical 
units; the actual displacements are always lower than 60 
h _1 kiloparsec, a value smaller than the typical core radius 
of the ICM distribution. 

A final test to see if the characteristics of the detector 
play a role in the offset measure can be finally performed 
on experimental data. If the measured offsets were pro- 
duced by some unknown instrumental effect, we would 
expect to find a correlation between the redshift of the 
cluster and the offset measured in physical units (kilopar- 
sec). The linear Pearson correlation coefficient r between 
the two dataset is r = 0.25. The probability that two un- 
correlated 22 elements datasets give a value of r > 0.25 is 
P = 40%. 



6. X— ray cD counterparts 

Once the presence of offsets has been firmly established, 
we can try to correlate both X-ray positions with the op- 
tical position of the dominant cluster galaxy. 

Unfortunately, the comparison between optical and X- 
ray positions suffer the boresight problem, i.e. an absolute 
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Fig. 5. Distribution of the measured offsets in physical units. 



error in the satellite pointing which turns in a system- 
atic shift between sources positions and respective optical 
counterparts. This error can be corrected when at least 2 
bright X-ray sources in a single image can be correlated 
with bright optical objects (e.g. Guide Star Catalog ob- 
jects). If uncorrected, the boresight error is estimated to 
have a Gaussian distribution with standard deviation of 
~ 10" both in right ascension and declination (see the 
ROSATSRC catalog, Zimmermann 1994). 

In the A478 image, after the application of a 12" bore- 
sight correction, the X-ray peak position is only 3" away 
from the cD optical position. 

No X-ray source correlated to a Guide Star Catalog 
object was found in the remaining 21 images and, given 
the rather large pointing error, it has not been possible to 
establish a firm identification of the cD galaxy as the ob- 
ject that produces the offset peak in the large scale X-ray 
emission of the ICM. However, we can state that within 
these uncertainties the position of the X-ray peak coin- 
cides with the galaxy. In figure |^ we show the distribution 
of X-ray peak positions around the cD which lies in the 
origin of the coordinates. The circles represent la and 2a 
confidence regions. As expected, 40% (9 over 22) of the 
X-ray peaks lie inside the la circle. The distribution of 
the Right Ascension residuals in figure |s| appears biased 
towards negative shifts. The tests performed on the anal- 
ysis procedure described in section |j shows that no bias is 
present in the X-ray positions and the accuracy of Guide 
Star Catalogue is fully adequate; moreover, the absence of 
the effect in the declination shifts indicates that the the 
analysis procedure is unbiased. The binomial probability 
of obtaining such a skewed dataset from a symmetric dis- 
tribution is P = 3.2%. 

If the difference between the peak and X-ray positions 
were not due to instrumental effects we would expect a 
correlation between angular offset and redshift. The corre- 
lation coefficient between these two datascts is r = —0.11. 
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Fig. 6. Relative positions of the X-ray peaks with respect to 
the cD galaxy. The cD is located in the origin of coordinates 
and shifts in right ascension and declination are given in arc- 
seconds. The dashed circle encloses the la confidence region 
for the boresight error while dotted circle marks the 2a region. 



The possible identification of the cD galaxy with the 
peak of the cluster emission has been discussed in several 
papers. Allen et al. (1995) find a coincidence of the two 
positions in clusters that host a massive cooling flow in 
their centre, but do not find evidence of centroid shift 
in their sample; more often, a coincidence between the 
orientation of the cD isophotes with the ICM emission 
major axis is found (Sarazin & McNamara 1997; Huang 
& Sarazin 1998). 

Given the identification, we can discuss a model that 
predicts the amount of separation between the X-ray peak 
and the bottom of the potential well of the ICM, marked 
by the centroid of large scale emission. 

7. The oscillatory model 

An implication of the above discussion is that the cD 
galaxy is not located in the bottom of the potential well 
and hence cannot be at rest but must oscillate around the 
equilibrium position. This wobbling produces the observed 
offset. 

With a small set of reasonable assumptions and/or ap- 
proximations on cluster structure, we can obtain the equa- 
tion of motion of the cD in the bottom of the potential 
well. First, since dominant galaxies are found well inside 
the core radius of clusters (section |J) , we assume that the 
mass inside the area of oscillation is uniformly distributed. 
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In clusters that host a massive cooling flow in their 
centres great deviations from uniformity could be present 
in the baryon fraction, however they will have a small 
effect on the total mass distribution and hence do not 
invalidate our model. We have therefore: 
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m,QX 



= F = -G 



Mi nt ma 



= --nG ptotm G x 



(3) 



where mo is the mass of the galaxy, Mi„ t the mass con- 
tained in a sphere of radius equal to the position of the cD 
and ptot the total mean mass density in the central part 
of the cluster. Solving equation 3 we obtain: 



x(t) = B cos(uit + p) 



where: 



u = \l ^ nG Ptot 



(4) 



(5) 



The actual values of Bj and ipj for the j-th cD galaxy 
are linked to initial speed and position. Since we are con- 
cerned with quantities averaged on an ensemble of clus- 
ters, we expect the values of the phases <p to be erased. The 
value of B can be obtained by assuming the rms speed of 
the galaxy proportional to the radial velocity dispersion 
of the whole cluster. This seems to be a reasonable as- 
sumption in view of the fact that the kinematic of the cD, 
whatever it is, must be related to the cluster potential 
well. 

Putting this relation in our model we have: 
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Fig. 7. Dependence of the physical offset from the radial 
velocity dispersion. The offset histogram of the high disper- 
sion subsample (< a v >= 888km/ s, dashed line and ar- 
row) is significantly larger than that of the low dispersion 
(< <J V >= 597fcm/s, solid line and arrow) one. The arrows 
mark the mean value of the two distributions. 



To test the likelihood of the assumption that the cD 
oscillation velocity is linked to the whole cluster velocity 
dispersion, made in equation 6, we have divided the sam- 
ple in two subsamples based on their radial velocity dis- 
persion. We expect the average offset to be proportional 
to the average velocity dispersion in the two subsamples. 
As is shown in figure |?], denoting with ^ and ; respectively 
the high and low velocity dispersion subsamples we have: 



= < v 2 > - < v > 2 =< v 2 > 
= B 2 u) 2 < sin 2 (cot) > 



from which: 



B 2 



2a 2 



2 nGpto 



(6) 



(7) 



We can hence define an adimensional oscillation am- 
plitude r normalized to the individual cluster properties 
that should be distributed as an ensemble of identical har- 
monic oscillators with the same amplitude but with ran- 
dom phases and seen from random orientations. 



r = 



p(T) cx 



2ttG czd$ p\gt 
3 Hq a v 



dl 



It ^(l 2 ~T*){V -P) 
where L is the largest oscillation amplitude allowed. 



(8) 



(9) 



< a v >h= 888 km/ s 

< a v >i= 597 km/s 

from which: 

< &v >h „„ „ < a v > 



< A > h 



37.9 



< A> 



< A > h = 2iAh~ 1 kpc 

< A > ; = W.lh^kpc 



37.1 



(10) 



The major problem in the application of equation^ is 
the determination of the central mean mass density. In 
fact, the X-ray emission of the ICM gives a direct infor- 
mation of the baryon density while what we need is an 
estimate of the total mean density in the cluster centres. 

A crude estimate of the mean central density can be 
inferred applying the virial theorem with the measured 
cluster parameters a v and r c : 



1 el 



1 ol 



Ptot ~ Q R2 ~ Q ^ 



r = 



Inserting this density in equation^ we obtain: 
27T czdtf 1 2ir d$ 



3 Hi 



(I 



3 r, 



(11) 



(12) 
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Fig. 8. Histogram of the values of the F parameter, as de- 
scribed in the text. The shaded area represents the expected 
distribution binned with the same width of data. A Kol- 
mogorov-Smirnov test gives a ~ 75% probability of success. 
The insert shows the theoretical distribution (eq. H) before the 
binning process. Poissonian errors are drawn on experimental 
points. 



where r c $ is the apparent angular core radius. 

A further problem that affect the comparison between 
the harmonic model and the observations is the intrinsic 
difficulty in measuring the core radius of clusters that host 
a cooling flow in their centre. This measurement is made 
highly inaccurate by the need of masking the central re- 
gion of the cluster to get rid of the modification on the 
cluster central surface brightness due to the emission of 
the cool and dense inflowing gas (see, e.g. Allen & Fabian 
1997). To see which uncertainty this problem could con- 
tribute to the final parameter T we have tested the change 
of the fitted value of the core-radius increasing the radius 
of the masked central region in the image of A1795. When 
the masked region is smaller than the PSF FWHM, the 
core radius is highly underestimated. If the mask radius is 
bigger than twice the PSF FWHM - any influence of the 
central compact source being completely erased - the core 
radius changes by a factor of 30% increasing the mask 
radius up to 5 times the above value. This uncertainty 
propagates with unchanged percentage in the determina- 
tion of the value of T. Given the smooth shape of the 
probability distribution (eq. ||) we expect that these un- 
certainties will affect only the high amplitude sharp cutoff 
and the low amplitude cusp of the function, leaving the 
overall shape unmodified (see the insert in figure ||). 

Figure |§| shows the distribution of the measured T in 
our cluster sample, compared with the theoretical distri- 
bution obtained rebinning equation^ (shaded area). The 
best value for the largest amplitude L = 0.69 has been ob- 



tained maximizing the Kolmogorov-Smirnov probability 
that the measured distribution is drawn from that theoret- 
ically derived. The maximum probability is Pk-s — 73%. 
Figure || may suggest a slightly bigger value for the max- 
imum amplitude L, however this would require a higher 
number of small offset clusters. This discrepancy can be 
ascribed to the crudeness of the model and, in particular, 
the introduction of circular (or even elliptical) cD orbits 
would erase the cusp of the distribution in the lowest off- 
set region. Such a refinement of the model would however 
require a larger dataset to allow any quantitative compar- 
ison. 



8. Conclusions and discussion 

We have analyzed the structure of the X-ray emission of 
a set of 22 Abell cD clusters. We find evidence for a dis- 
placement between the centroid of the large scales and the 
peak of the emission, the latter being coincident with the 
dominant galaxy position. 

The presence of centroid shifts in clusters with cD 
galaxies and cooling flows is problematic since such off- 
sets are expected to be due to subcluster mergers, phe- 
nomena that would disrupt cooling flows (Meiksin 1988; 
Friaca 1993). 

We have shown that, identifying the peak with the cD 
galaxy, an harmonic oscillator model for its motion in the 
cluster potential well is in agreement with present obser- 
vations and with the predictions of the numerical simula- 
tions (Malumuth & Richstone 1984). It follows that offset 
peaks in X-ray images are produced by this oscillation 
and should not be thought as signs of dynamical instabil- 
ity and/or substructure. 

From this interpretation, it follows that the compact 
emission in the centre (or near the centre) of clusters can 
be produced in principle inside the central galaxy or by 
a fraction of the ICM strongly linked to the galaxy itself. 
The first hypothesis can be easily discarded since obser- 
vations of cluster central regions with high resolution in- 
struments show that the compact source has a small but 
non-zero extension and hence cannot be due to an active 
nucleus (see, e.g., Grebenev et al. 1995, Pierre & Starck 
1998). In the second scenario the cool and dense gas of 
the cooling flow could be responsible for the peaked emis- 
sion seen in X-ray images and the offset would be caused 
by the motion of the cD galaxy. In this case the cooling 
flow - or at least its inner part - would have to follow the 
cD galaxy in its motion. Whether this is possible or such 
a motion would disrupt the frail flow stability has to be 
investigated in more detail. From the observational point 
of view, Pierre & Starck (1998) have analyzed a sample of 
high resolution images of clusters, finding that in the inner 
cores of massive cooling flow clusters the X-ray emission 
shows peculiar features and strong isophotal twisting (see 
their figure 6). 
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As a further consequence of this interpretation, the 
fraction of structured clusters is greatly reduced. In fact, 
in his analysis of a flux limited sample of X-ray clusters, 
Davis (1994) found a residual emission near the centre of 
clusters (but not coincident with them) after subtraction 
of the extended component with the elliptic isophotal fit- 
ting technique. The interpretation of this feature as a sign 
of substructure led to the conclusion that up to 70% of the 
X-ray clusters were structured. If their results are reinter- 
preted on the light of the conclusions of this work, we find 
that the above fraction falls to less than 30%. 

The hardest problem we had to deal with in this work 
is represented by pointing errors in X-ray images. The 
boresight uncertainties do not allow us to place tight con- 
straints on any significant displacement between the peaks 
of the X-ray and optical emission and, even if figure |^ im- 
plies the positions are in agreement, any improvement will 
require better astrometric accuracy through new observa- 
tions. Note, however, that the comparison with numerical 
simulations and the calculation of the parameter T de- 
scribed in section |?] have been performed only using rel- 
ative positions from the X-ray images and hence do not 
suffer of any systematic pointing error. 

The new generation of X-ray telescopes (JET-X, 
AXAF, XMM, WFXT, ...), thanks to their higher spa- 
tial resolution and sensitivity, will presumably give us the 
possibility of better testing the model. Moreover, it is de- 
sirable that new observations will allow to enlarge the set 
of cluster observations to be compared to the predicted 
harmonic amplitude. 
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